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摘 要 : 操作 工人 在 长 时 间 的 工作 过 程 中 ， 由 于 精神 和 身体 的 疲劳 会 产生 失误 ， 人 的 可 靠 性 降低 ， 导 致 生产 残 次 品 ， 
甚至 危及 人 身 安全 。 为 了 了 解 在 保证 可 靠 性 的 前 提 下 ,操作 者 最 长 可 持续 的 工作 时 间 ( 极 限 工作 时 间 )， 提 出 了 一 种 结合 
统计 过 程控 制图 (statistical process control chart，SPCC) 和 差分 自 回归 移动 平均 模型 (autoregressive integrated moving 
averagemodel，ARIMAI) 的 方法 来 进行 确定 。 用 连续 时 间 中 单位 时 间 间 隔 内 的 次 品 率 衡 量 人 的 可 靠 度 ， 并 形成 次 品 率 时 
序 序列 (defective rate time series, DRTS)。 一 方面 使 用 该 数列 基于 ARIMA 构建 SPCC, 并 确定 样本 的 上 控制 线 (up control 
line, UCD)、 下 控制 线 (low controlline,LCL) 和 极限 时 间 范 围 ; 另 一 方面 使 用 ARIMA 和 支持 向 量 机 (support vector machine， 
SVM) 配 合 实时 跟踪 算法 (realtime tracing algorithm，RTA) 对 需要 确定 极限 工作 时 间 (limit working time，LWT) 的 某 一 操 
作者 的 DRTS 进行 预测 。 最 后 将 预测 曲线 置 于 SPCC 中 来 确定 该 操作 者 LWT 的 范围 。 通 过 对 10 个 操作 者 8 小 时 ， 每 
陪 10 分 钟 采 样 得 到 的 DRTS 分 析 ， 表 明 该 类 工作 参数 的 范围 为 UCL=[3.89，4.66]、LCL=[-4.98，-3.90]、LWT=[393， 
450] min。 确 定 了 一 个 新 操作 者 的 LWT 为 [435，452] min， 即 435 min 时 该 操作 者 就 应 该 停止 工作 进行 休息 。 
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Abstract: In the process of long work hours, operators can produce errors due to mental and physical exhaustion, human 
reliability is reduced, and it leads to production of incomplete and even endangers humanal safety. In order to understand 
operator's longest sustainable working time on the premise of given reliability, a method of combining SPCC and ARIMA is put 
forward. With unit time interval defect rate in the continuous time to measure the reliability of human, and form a defective rate 
time series. On the one hand, the SPCC is built using DRTS based on ARIMA to determine the sample of UCL, LCL and LWT 
range. On the other hand, the DRTS of an operator is predicted using ARIMA and SVM with RTA. Finally the prediction curve 
is put into the SPCC to determine the LWT range. Based on 8 hours, 10 operators DRTS analysis of sampling for every 10 
minutes, indicates that the scope of this kind of working parameters for the UCL = [3.89, 4.66], LCL = [4.98, 3.90], LWT = 
[393450]. Identified a new operator LWT is [435452] minutes, 435 minutes when the operator should stop working to have a 
rest. 


Key words: reliability of the human; ARIMA; SPCC; limit Working time 


心情 烦躁 ， 注 意 力 难以 集中 ， 不 但 工作 效率 低下 ， 而 且 会 生产 

加 较 多 残 次 品 。 这 样 会 加 重 操作 人 员 的 心里 压力 ， 使 人 的 行为 更 
我 国 是 一 个 劳动 力 大 国 ， 在 高 度 机 械 化 和 自动 化 的 今天 ， 加 不 可 靠 ， 甚 至 会 由 于 操作 不 当 造 成 人 身 伤 亡 。 所 以 ， 确 定 一 

仍 有 大 量 的 劳动 者 从 事 着 以 人 为 主要 动力 来 源 的 劳动 。 人 与 机 ”个 适当 的 工作 时 间 ， 不 仅 能 提高 操作 者 的 效率 ， 而 且 能 提高 产 
器 相 比 有 着 丰富 的 创造 能 力 ， 但 同时 人 在 长 时 间 的 劳动 过 程 中 品 的 合格 率 。 一 方面 保护 了 劳动 者 ， 另 一 方面 也 增加 了 企业 的 
会 在 生理 和 心里 方面 产生 疲劳 ,在 疲劳 状态 下 工作 的 操作 人 员 ， ”经 济 效 益 。 即 在 保证 操作 人 员 操 作 可 靠 性 不 小 于 某 个 特定 值 的 
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前 提 下 ， 确 定 操作 者 可 以 工作 的 最 长 时 
(LWT)。 

关于 人 的 可 靠 性 及 适合 的 工作 时 间 的 研究 ， 目 前 有 柴 松 等 
人 [2 的 基于 CREAM 和 不 确定 推理 的 人 因 可 靠 性 分 析 方法 和 
海洋 工程 人 因 可 靠 性 定量 分 析 方法 与 应 用 ， 董 学 军 等 四 的 基于 


司 ， 即 极限 工作 时 间 
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洽 铁 军 ， 等 : 生产 过 的 极限 工作 时 间 研 究 
从 图 1 可 知 算法 分 为 左右 两 部 分 用 粗 虚线 分 开 。 左 边 为 时 


序数 据 的 预测 部 分 ， 右 面 为 SPCC 图 构造 部 分 。 最 右 侧 的 
ARIMA 处 理 模 块 包括 了 三 个 部 分 ， 左 右 两 部 分 都 要 用 
独立 开 来 , 右 侧 形成 SPCC 图 , 是 来 源 于 多 个 样本 的 统计 结果 ， 
是 整个 算法 的 基础 。SPCC 可 得 到 LCL、UCL、CL 和 LWT， 


补偿 和 不 可 蔡 代 因素 合成 的 人 因 可 靠 性 分 析 方法 ， 蒋 建 军 等 人 
5] 的 基于 隐 马 尔 可 夫 的 核电 厂 半 数字 化 人 -机 界面 事故 诊断 过 
程 人 因 可 靠 性 模型 和 考虑 人 因 的 核电 厂 主 控 室 认 知 可 靠 性 模型 
研究 ， 侦 鞭 外 的 煤矿 职工 疲劳 度 相 关 因素 分 析 与 控制 研究 。 

笔者 参考 了 文献 [7~9]， 并 根据 自己 的 经 验 00I9， 提 出 了 基 
于 SPCC 和 ARIMA 的 人 的 LWT 的 确定 方法 。 该 方法 主要 分 为 
两 部 分 , 一 部 分 是 构造 SPCC 并 得 到 LCL、UCL、 中 心 线 (center 
line，CL) 和 极限 时 间 范 围 ， 这 部 分 的 数据 来 源 是 多 个 样本 ( 操 
作者 ) 的 统计 数据 ， 男 一 部 分 使 用 ARIMA， 配 合 SVM 和 RTA 
预测 某 一 个 新 样本 在 工作 时 间 不 足 的 情况 下 后 继 时 间 序 列 的 变 
化 ， 进 而 形成 该 新 样本 的 DRTS， 结 合 得 到 的 SPCC 可 预测 该 
样本 的 LWT。 
1 ”方法论 

该 方法 涉及 到 的 相关 知识 有 SPCC、ARIMA,SVM 和 RTA。 
算法 作用 有 两 个 ， 一 是 统计 多 个 样本 操作 者 〉 确定 某 种 操作 
的 泛 化 LCL、UCL、CL 和 LWT; 二 是 基于 这 些 统计 结果 预测 
某 个 个 体 ( 操 作者 ) 的 LWT。SPCC 对 多 个 样本 统计 数据 进行 
整理 ， 得 到 数据 的 LCL、UCL、CL 和 LWT。 可 得 数据 变化 趋 
势 是 否 将 超 限 ， 从 而 可 根据 统计 数据 的 特征 ， 预 测 同类 型 数据 
是 否 有 超 限 的 可 能 , 形成 SPCC 图 是 衡量 标尺 。 ARIMA 是 算法 
核心 ， SPCC 处 理 时 间 序 列 。ARIMA 在 形成 SPCC 时 适用 于 无 
相关 性 的 时 间 序 列 数据 。 为 了 预测 新 操作 者 的 LWT， 需要 对 短 
时 间 内 的 新 操作 者 的 DRTS 进行 分 析 ， 并 对 未 来 的 序列 数据 进 
行 预测 。 但 DRTS 很 可 能 是 线性 和 非 线 性 的 序列 数据 县 加 的 ， 
根据 文献 [9] 可 将 线性 数据 使 用 ARIMA 进行 预测 ， 非 线性 数据 
使 用 SVM 预测 .根据 文献 [8] 可 使 用 RTA 对 时 序数 据 进 行 跟踪 ， 
以 保证 预测 的 精确 性 。 经 过 上 述 处 理 的 线性 和 非 线 性 时 序数 据 
相 县 加 形成 最 终 的 该 操作 者 的 DRTS 数据 ， 放 入 SPCC 图 中 ， 
即 可 确定 他 的 LWT。 算 法 流程 如 图 1 所 示 。 
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图 1 算法 的 流程 


进而 处 理 单个 样本 的 LWT 预测 。 左 侧 为 预测 单个 样本 DRTS 数 
据 ， 在 SPCC 中 与 LCL 和 UCL 比较， 最终 确定 LWT。 


2 ”方法 步骤 及 其 解释 


SPCC 对 实时 数据 进行 监控 和 预测 。SPCC 基于 ARIMA 与 
传统 的 SPCC 是 不 同 的 0544。SPCC 构造 过 程 分 三 步 : 相关 性 
平价 ; ARIMA 构建 ; 控制 图 构建 。 自 相关 函数 (autocorrelation 
coefficient function，ACF) 在 第 一 步 中 ， 用 来 评价 时 间 序 列 的 自 
相关 性 。 如 果 不 相 关 直 接 执行 第 三 步 ， 用 原始 时 间 序 列 构 建 控 
制图 。 和 否则 用 ARIMA 消除 相关 性 ， 在 第 二 步 产 生 时 间 序列 的 
非 相关 残 差 。 第 三 步 用 时 间 序 列 残 差 构建 控制 图 。 
2.1 相关 性 估计 
使 用 ACEF 估计 相关 性 。 设 滞后 时 间 间 隔 为 k， 使 用 
MATLAB 中 的 autocorr 计算 相关 性 。 如 果 ACEF 图 显示 的 数值 
绝对 值 总 体 趋 势 减 小 ， 则 断定 时 间 序 列 是 相关 1。 
2.2 ARIMA 构建 

如 果 数 据 是 相关 的 , 那么 使 用 ARIMA 消除 来 消除 相关 性 。 
ARIMA 模型 可 表示 为 ARIMA(p,d,q)， 其 中 p,d,q 为 大 于 等 于 0 
的 整数 , p 为 自 回 归 项 ; d 为 时 间 序 列 成 为 平稳 时 所 做 的 差分 次 
数 ，q 为 移动 平均 项 数 。 给 定 ARIMA(p,d,q)， 原 时 间 序 列 x 可 
被 拟 合并 转换 为 一 个 白 噪 声 过 程 @。, 之 间 的 关系 如 式 (1) 所 示 
[07] 


< 


-S08 as -0+ Spa 上 
p=1 


a=l 


这 


其 中 ，W 和 gp, 自 回归 系数 和 移动 平均 系数 ，B 延迟 算 子 ， 
BX = 7 为 延迟 次 数 。 
ARIMA 的 构造 分 三 个 步骤 , 即 识别 规则 , 参数 估计 和 模型 
检验 03。 
a) 在 识别 规则 中 ， 使 用 逐次 差分 将 时 间 序 列 转换 为 稳定 数 
据 ， 一 次 和 二 次 差分 公式 分 别 如 式 (2)(3) 所 示 趾 。 
(1—B)x=%—X (2) 
(1-B)%=(-2B+B)x, =% 2X +N; G3) 
ACF 和 偏 自 相关 函数 (partial autocorrelation coefficient 
fonction，PACF) 可 绘制 出 各 种 差分 时 间 序 列 的 对 应 曲线 。 如 果 
其 他 参数 不 变 , 差分 次 数 d 的 变化 对 ACF 和 PACF 形成 的 图 
影响 较 小 091]。 有 具体 来 说 ，ACF 图 像 拖 尾 和 PACF 延迟 p 后 截 
时 应 使 用 ARIMA， 具 体 方法 见 文献 [19]。 
b) 使 用 最 大 化 条 件 对 数 似 然 函 数 7(g,,0,,2) 进行 估计 ， 如 
式 (4) 所 示 中 。 


是 式 


S:(9,,0,) 
202 


14200550)=-5In() (4) 
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其 中 ， 是 e 的 方差 ，S.(p,,@ ) 通过 式 (5) 进 行 估计 。 


S$.(9,,0,) = De:(p,,0, |0 -Bx ,e,,(1—B)x, ) (5) 


如 果 ，ARIMAC,dq) 中 pz 上 0,g# 上 0， 通 过 最 小 艾 卡 信息 量 
准则 (lowest Akaike information criterion，AIO) 进 行 定 阶 20， 其 
具体 表达 式 如 式 (6) 所 示 。 


AIC(, )) =In(6? )+ 了 (0) 


其 中 : 5。 是 6 的 方差 的 估计 值 ，N 为 样本 数 。 


0o) 检 验 时 间 序 列 的 相关 性 是 否 消除 。 使 用 ACF 对 白 噪 声 e 
进行 分 析 可 进行 验证 。 当 延迟 7 大 于 10, 且 ACF 的 分 析 结 果 趋 
近 于 0 时 可 认为 相关 性 已 被 消除 PCJ。 和 否则 重复 进行 ARIMA 处 
里 ， 直 到 满足 上 述 条 件 。 
d) 利 用 RTA 改进 ARIMAEI。 在 ARIMA 基础 上 提出 与 实时 
跟踪 算法 RTA 结合 的 改进 模型 。 假 设 要 对 时 间 序 列 {x |i=1:7} 
进行 预测 ， 假设 步 长 为 m， 利 用 {x|i=1:t} 预测 
{w|i=t+1:t+m} ， 去 掉 {x|i=1:1+m} 的 前 m 个 数据 ， 形 成 
{x | 上 =m+1:t+m} 进行 下 一 步 的 时 序数 列 预测 。 

e) ARIMA 和 SVM 模型 分 别 在 处 理 线性 和 非 线 性 数据 结构 
的 预测 问题 时 有 着 各 自 的 优势 。 可 将 ARIMA 和 SVM 结合 使 
用 ， 以 求 达 到 增强 预测 模型 对 数据 样本 的 适应 能 力 ， 从 而 提高 
预测 的 效果 外。 

2.3 控制 图 构建 

构建 体 哈 特 统计 控制 过 程 图 (Shewhart statistical process 
control chart, SSPCC), 基于 Shewhart three sigma control theoryP5] 
来 确定 中 心 线 CL。 历 史 数据 可 为 无 相关 性 的 原 时 间 序 列 x 或 


ARIMA 处 理 后 的 白 噪声 过 程 。。@& 为 在 时 间 点 i 处 的 历史 数 
据 。 控 制 上 线 UCL、 控 制 下 线 LCL 和 中 心 线 CL 如 式 (7) 所 
示 。 


UCL/ LCL=2CL =X+30=X+3™ (7) 
其 中 。 。 宛 ，T 为 时 间 序列 长 度 的 平均 值 。 二 表示 平均 
-与 - 
移动 范围 ， 全 Wi ， 其 中 jm =@ 一 ai。 
mr = 三 
7-1 
SSPCC 使 用 时 间 序 列 中 两 个 连续 值 作为 移动 范围 2, 所 以 


n=2, qd, =1.128 1。 
3 ”实例 计算 


在 流水 线 上 的 操作 人 员 往 往 由 于 工作 时 间 长 而 产生 疲劳 ， 
导致 生产 的 部 件 残 次 品 增多 ， 甚 至 可 能 由 于 操作 不 当 对 自身 安 
全 遭 威 胁 。 监测 10 个 工人 在 8 小 时 内 的 生产 残 次 品 比例 (次 品 
率 ) 作为 衡量 操作 人 员 操 作 可 靠 性 指标 。 取样 时 间 间 隔 为 10 分 
钟 ， 次 品 率 为 10 分 钟 内 生产 的 残 次 零件 /全 部 零件 。 构 成 了 
48 个 数据 点 的 10 个 样本 时 序 序列 数据 。 根 据 次 品 率 统计 ， 在 
8 小 时 内 的 10 个 工人 的 次 品 率 在 0~0.2， 曲 线 是 起 伏 的 ， 但 总 
趋势 递增 。 所 以 采用 一 阶 差 分 对 时 序 序 列 进行 处 理 。 由 于 一 阶 
差分 、.ACF、PACEF 在 MATLAB 中 已 有 工具 包 , 这 里 不 在 说 明 。 
给 出 10 个 样本 中 第 一 个 样本 的 数据 处 理 过 程 的 图 示 , 如 图 2 所 
不 。 
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2 中 , (a) 表 示 了 该 操作 者 (第 一 个 样本 ) 的 DRTS 数据 ， 


图 2 第 一 个 样本 的 数据 处 理 过 程 图 


图 (和 (可 见 ， 大 部 分 数据 的 ACF 和 PACF 值 都 超出 允许 范 


一 共 48 个 采样 点 ，480 分 钟 的 工作 时 间 。 根 据 次 品 率 定义 ， 本 
知 DRTS 数据 是 起 伏 波动 的 ， 随 时 间 增 加 总 体 增加 的 趋势 。 


姑 ， 时 间 序 列 的 平稳 性 不 满足 计算 要 求 ， 需 要 进行 差分 计算 。 
通过 差分 计算 的 一 次 差分 序列 .ACF 和 PACF 分 别 如 图 (b)(d) 和 
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(e) 可 见 一 次 差分 序列 值 趋 于 稳定 ，ARIMA 模型 中 的 系数 4 定 
为 1。 然 后 利用 AIC 来 确定 ARIMA 模型 的 阶 数 p 和 q。 
计算 发 现 当 p=4,q=4 时 AIC 值 最 小 , 故 使 用 模型 ARIMA(4,1,4)。 
第 一 个 操作 工人 的 差分 DRTS 的 SPCC 如 图 3 所 示 。 


统计 过 程控 制图 
ee eo ee 
8. 
一 4 UCL 
让 
好 2r 
是 CL 
AR 0 
站 
-2 
-4 | | » | i | JILCL 
02468101214161820222426283032343638404244464850 
取样 点 /10 分 钟 


图 3 第 一 个 操作 工人 的 差分 次 品 率 时 序 序列 的 SPCC 


从 图 3 可 知 第 一 个 操作 工人 的 CL=0、UCL=4.17、LCL=- 
4.17。 差分 DRTS 曲线 与 UCL 交 于 采样 点 43~44, 即 在 430~440 
min 该 操作 人 员 达 到 限 值 , 不 适宜 继续 工作 。 则 LWT 范围 就 是 
430~440 min (LWL=437 min )。 同 理 确 定 10 个 操作 人 员 的 LWT 
如 表 1 所 示 。 将 表 1 中 参数 的 值 进行 综合 得 到 了 该 类 工作 的 对 
应 参数 的 范围 ， 即 UCL=[3.89,4.66] 、LCL=[-4.98,-3.90] 、 
LWT=[393,450] min。 


表 1 10 个 操作 人 员 的 极限 工作 时 间 


指标 1 2 3 4 5 6 7 8 9 10 范围 


UCL 4.17 4.08 3.89 5.11 4.66 4.54 3.90 4.43 3.68 3.99 [3.89,4.66] 


LCL -4.17 -4.08 -4.03 -4.98 -4.66 -4.34 -3.90 -4.33 -4.09 -4.18 [-4.98,-3.90] 


CL 0 0 -0.07 0.065 0 0.1 0 0.05 -0.205 -0.095 


LWT/ 分 437 404 401 450 422 430 393 425 397 406 [393,450] 


注 :范围 确定 是 取 10 个 样本 对 应 值 的 最 大 和 最 小 值 ,作为 区 间 的 两 个 端 


io 


上 述 过 程 完成 了 图 1 中 右 侧 框图 部 分 ， 即 完成 了 SPCC 的 
构建 。 下 面 是 判断 某 一 新 样本 《新 操作 者 ) 的 LWT。 
对 于 一 个 新 操作 者 〈 不 在 样本 的 10 个 操作 者 内 )， 利 用 以 
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为 2.459%)。 而 且 使 用 预测 曲线 结合 以 得 到 的 UCL 和 LCL 区 
域 可 确定 该 操作 者 的 LWT。 如 图 中 的 坚 向 底 纹 区 域 , 该 区 域 的 
水 平方 向 区 间 宽 度 是 由 DRTS 的 预测 曲线 与 UCL 相交 的 区 域 
宽度 确定 的 。 对 于 该 操作 者 的 LWT 为 [435,452] 分 钟 , 即 435 分 
钟 时 该 操作 者 就 应 该 停止 工作 进行 休息 。 
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图 4 次 品 率 时 序 序列 数据 差分 拟 合 的 残 差 


4 UCL= 


| [3.89.4.66] 
LCL= 
[-4.98,-3.90] 
取样 ng ) 钟 
图 5 该 操作 者 的 LWT 分 析 图 
4 ”结束 语 
a) 提出 了 极限 工作 时 间 确 定 方法 。 分 为 两 部 分 , 一 是 构造 


att 


SPCC 并 得 到 LCL、UCL、CL 和 极限 时 间 范 围 ， 数 据 来 源 是 多 
个 样本 (操作 者 ) 的 统计 ; 二 是 使 用 ARIMA, 配合 SVM 和 RTA 
预测 某 一 个 新 样本 在 工作 时 间 不 足 的 情况 下 后 继 时间 序 列 的 变 
化 ， 进 而 形成 该 新 样本 的 DRTS。 结 合 得 到 的 SPCC 可 预测 该 
样本 的 LWT。 
b) 基于 10 个 样本 形成 了 SPCC。 使 用 模型 ARIMA(4,1,4) 
进行 SPCC 构建 。 得 到 了 10 个 样本 (操作 者 ) 的 UCL、LCL、 


进行 工作 时 间 的 DRTS 的 一 次 差分 序列 对 该 操作 者 进行 LWT 
预测 。 利 用 文献 &9 的 方法 进行 时 序 序列 的 拟 合 和 预测 ， 得 到 了 
较 好 的 预测 结果 。 主 要 思想 是 将 原始 时 序 序列 拆 分 为 线性 序列 
和 非 线性 序列 ， 线 性 序列 由 ARIMA 处 理 并 进行 预测 ， 处 理 过 
程 与 上 述 过 程 相近 , 见 参考 文献 由 ; 非 线 性 序列 可 由 SVM 处 理 
进行 预测 ， 见 参考 文献 外。 同时 在 处 理 过 程 中 结合 RTA 方法 
进行 数据 跟踪 ， 提 高 预测 精度 。 最 后 将 上 述 线性 和 非 线性 序列 
进行 个 加 ， 达 到 对 原 DRTS 数据 预测 的 效果 。 实 际 值 与 预测 值 
的 残 差 如 图 4 所 示 。 预 测 误差 为 2.459% ,该 操作 者 的 实际 DRTS 
数据 曲线 和 预测 曲线 如 图 5 所 示 。 

在 图 5 中 , 底 纹 区 域 UCL 和 LCL 是 通 
定 的 。 使 


过 原始 10 个 样本 确 
预测 方法 对 该 操作 者 的 DRTS 预测 较为 成 功 〈 误 差 


Pm 


CL 和 LWT。 综 合 得 到 了 该 类 工作 的 对 应 参数 范围 ， 即 
UCL=[3.89,4.66]、LCL=[-4.98,-3.90]、LWT=[393,450] min。 

c) 确定 了 新 样本 (新 操作 者 ) 的 极限 工作 时 间 。 预 测 误差 
为 2.459%， 确 定 了 该 操作 者 的 LWT 为 [435,452] min， 即 435 
min 时 该 操作 者 就 应 该 停止 工作 进行 休息 。 
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